Astronomy & Astrophysics manuscript no. 0693text 


February 2, 2008 


(DOI: will be inserted by hand later) 





The role of the time step and overshooting in the modelling of 
PMS evolution: the case of EK Cephei 

J. P. Marques 1 - 2 , J. Fernandes 2 ' 3 , and M. J. P. F. G. Monteiro 1 ' 4 

1 Centra de Astrofisica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 

2 Grupo de Astrofisica da Universidade de Coimbra, Observatorio Astronomico da Universidade de Coimbra, Santa Clara, 
Coimbra, Portugal 

3 Departamento de Matematica da FCTUC, Coimbra, Portugal 

4 Departamento de Matematica Aplicada, Faculdade de Ciencias da Universidade do Porto, Portugal 
Received ; accepted 

Abstract. EK Cephei (HD 206821) is a unique candidate to test predictions based on stellar evolutionary models. It is a 
double-lined detached eclipsing binary system with accurate absolute dimensions available and a precise determination of the 
metallicity. Most importantly for our work, its low mass (1.12 M@) component appears to be in the pre-main sequence (PMS) 
phase. 

We have produced detailed evolutionary models of the binary EK Cep using the CESAM stellar evolution code (Morel fT997> . 
A ^-minimisation was performed to derive the most reliable set of modelling parameters (age, a A , a B and Y\). We have found 
that an evolutionary age of about 26.8 Myrs fits both components in the same isochrone. The positions of EK Cep A and B in 
the HR diagram are consistent (within the observational uncertainties) with our results. 

Our revised calibration shows clearly that EK Cep A is in the beginning of the main sequence, while EK Cep B is indeed a 
PMS star. Such a combination allows for a precise age determination of the binary, and provides a strict test of the modelling. 
In particular we have found that the definition of the time step in calculating the PMS evolution is crucial to reproduce the 
observations. A discussion of the optimal time step for calculating PMS evolution is presented. 

The fitting to the radii of both components is a more difficult task; although we managed to do it for EK Cep B, EK Cep A has 
a lower radius than our best models. 

We further studied the effect of the inclusion of a moderate convective overshooting; the calibration of the binary is not signifi- 
cantly altered, but the effect of the inclusion of overshooting can be dramatic in the approach to the main sequence of stars with 
masses high enough to burn hydrogen through the CNO cycle on the main sequence. 
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1. Introduction 

It is very well known that the interior structure of an evolved 
star is fixed by the initial mass M*, the initial chemical com- 
position (initial abundances in hydrogen, helium and metals, 
X[, Y{ and Z;, respectively) and the age f*. So, except for the 
Sun, modelling a single star using stellar evolutionary models 
on the HR diagram is not a closed problem because the number 
of parameters to be determined is larger than the observational 
constraints, the present observed luminosity and effective tem- 
perature. 

This difficulty is increased due to uncertainties on the 
physics of the models used to describe the stellar inte- 
rior/external structure. Currently these uncertainties are mod- 
elled by free parameters such as the diffusive and mass loss 
coefficient or the mixing length and overshooting parameters 
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for the convection in the MLT approximation (Mixing Length 
Theory). 

Some stellar multiple systems are in excellent position to 
avoid these difficulties. Good determinations of stellar masses 
are possible for the components of a given binary. Based on the 
reasonable assumption of a common origin for both compo- 
nents (which yields the same initial chemical composition and 
age), the problem of modelling both stars of a binary system 
is reduced to the determination of three so-called stellar mod- 
elling parameters, namely, Y u Z;, f* (since X+Y+Z-l we do not 
need to determine X[ independently). Moreover, if the stars are 
young enough (so that diffusion processes have not had time to 
change the surface composition significantly) and the present 
surface metallicity of the stars is also known (and is, therefore, 
equal to the initial metallicity), then (Z/X)i is known. So, the 
number of observables (five: the effective temperatures and lu- 
minosities for both stars plus the metallicity of the system) is 
higher than the number of modelling parameters. This allows 
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Table 1. Properties of EK Cep (HD 206821). All data except 
chemical composition are from Andersen ( I199U and Popper 
( 1987 »; chemical compositions are from Martin & Rebolo 
(FEU . See also Tomkin (TW51 . Ebbighausen fT%6l . and Hill 
& Ebbighausen (THJPl . 





EK Cep A 


EK Cep B 


M* /Mq 


2.029±0.023 


1.124±0.012 


log7/ eff ,* (K) 


3.954+0.010 


3.756+0.015 


log (LJLq) 


1.17±0.04 


0.21±0.06 


RJRe 


1.579±0.007 


1.315±0.006 


logg (cgs) 


4.349+0.010 


4.251 ±0.006 


Sp. Type 


A1.5V 


G5Vp 


[Fe/H] 




+0.07±0.05 



us to include unknown parameters related to the description 
of the interior physics. For low mass stars, without mass-loss 
or strong rotation, the convection is the main source of uncer- 
tainties in the description of the stellar interior. So we choose 
to include the mixing length parameter for each component of 
a binary system a A and a s . A number of binary systems have 
been calibrated in this way, e.g. a Cen (Noels et al.[T991 , Morel 
et al. l2000b> and t Peg (Morel et al.EDOOa); all these stars are, 
however, main sequence (MS) or post-main sequence stars. 

A good candidate binary system for calibration must have 
well-determined luminosities, effective temperatures, metallic - 
ities and dynamical masses. Pre-main sequence (PMS) binaries 
that have all the required characteristics are very rare at present 
(see Palla & Stahler 2001 for some examples of PMS binaries, 
and Lee et al. 1994 for lithium abundances for some of them). 

The double-lined eclipsing binary system EK Cep 
(HD 206821) has these unique characteristics. Accurate abso- 
lute dimensions are available (Andersen 1 199 II Popper |l 987 1 
and a determination of the surface metallicity was made 
(Martin & Rebolo 1993 1. Most importantly, the radii of the two 
stars are much closer than expected given their mass ratio. This 
indicates that EK Cep B is still contracting towards the main- 
sequence. At the same time, the higher mass of the primary 
implies that the star must be already in its main sequence evo- 
lution. All these make EK Cep a perfect candidate to test theo- 
retical models of pre-ZAMS solar-type stars. A summary of the 
observational characteristics of this binary is given in Tabled 

Several authors have tried to model this system. A prob- 
lem frequently found is that EK Cep A has a smaller radius 
(or that EK Cep B has a bigger radius) than that predicted by 
the models, making it difficult to fit both components to the 
observations in the same isochrone. Claret et al. (1995 1 have 
computed models that agree reasonably well with the obser- 
vations; however, the radius of EK Cep B is somewhat higher 
than their models predict (a problem also found in the models 
used for comparison with the observations in Martin & Rebolo 
1993 1. Yildiz (2003 1 computed models using a higher metal- 
licity than Claret et al. and a fast rotating core for EK Cep A 
in order to lower its luminosity and radius; this way, models 
that fit the observed radii and luminosities of the EK Cep sys- 



tem can be computed. The metallicity used in Yildiz (2003 
Z ~ 0.04) is much higher, however, than that given in Martin 
& Rebolo ( 119931 Z ~ 0.02). The values for the iron abundance 
given by these authors spread somewhat, making a very precise 
determination of this abundance difficult; nothing in these val- 
ues comes close, however, to the value of metallicity used in 
Yildiz J20031 

Our atempt to model this system focuses, instead, on the 
physical input that might affect the PMS evolution. Our own 
PMS models, computed using the CESAM stellar evolution 
code (Morel 1997 1, suggested that there are effects in the PMS 
evolution due to the physical inputs (overshooting, time step) 
that have been ignored or too quickly dismissed in the previous 
works. 

2. Physical ingredients 

Models were computed using the CESAM stellar evolution 
code (Morel 1997 1. Typically, each model is described by about 
600 shells, and an evolution by about 400-600 models. We re- 
stricted the maximum time step of the evolution: for EK Cep A, 
the maximum time step was kept at about 0.017 Myrs, while for 
EK Cep B we allowed the time step to reach 0.075 Myrs. We 
justify these values for the time step in section[6] 

Each evolution is initialised with an homogeneous, fully 
convective model in quasi-static contraction (Iben 1965 1, with 
a central temperature inferior to the ignition temperature of 
the deuterium. We shall call the "age" of the model the time 
elapsed since initialisation. This initial model for PMS evolu- 
tion is not a realistic assumption, since stars do not form by 
the homologous contraction of the pre-stellar cloud. Instead, 
a hydrostatic core first forms, which accretes mass from the 
parental cloud until either stellar winds disperse the cloud or 
the cloud material exhausts. Then the star becomes optically 
visible for the first time, since no more surrounding material 
obscures it: the star is born (see Stahler 1983, Palla & Stahler 
119911 I1992> . If we wanted to model very young stars, the 
problem of the initial conditions would be very important; for 
EK Cep, however, the initial conditions are irrelevant since it is 
much older than 1 Myr. Our "age" of the model, although not 
a real age of the star, is very close to it. 

The model of the zero-age main-sequence (ZAMS) is de- 
fined as the first model where nuclear reactions account for 
more than 99% of the energy generation. 

We used the OPAL equation of state (Rogers et al. 1996; 
we used the 2001 version of the OPAL tables which goes down 
to 2000 K) and the opacities of Iglesias & Rogers (19961, com- 
plemented, at low temperatures, by the Alexander & Fergusson 
dl994l l opacities. 

The temperature gradient in convection zones is computed 
using the standard mixing-length theory (Bohm-Vitense 1958 1. 
The mixing length is defined as l=a*Hp, Hp being the local 
pressure scale height, Hp- - dr/dlnP. 

The nuclear network we use contains the following species: 
'H, 2 H, 3 He, 4 He, 7 Li, 7 Be, I2 C, 13 C, 14 N, 15 N, lfi O, n O, 9 Be, 
and an "extra" fictitious non-CNO heavy element which com- 
plements the mixture; this element has atomic mass 28 and 
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charge 13. Deuterium and lithium burning are taken into ac- 
count, as well as the most important reactions of the PP+CNO 
cycles. The nuclear reaction rates are taken from the NACRE 
compilation (Angulo et al. 1999). 

Given the expected young age of the binary, we do not 
consider diffusion processes. However, these effects must be 
taken into account if we would follow the evolution of the sur- 
face chemical composition. If the surface chemical composi- 
tion were known with high accuracy for both components, the 
binary would be an excellent candidate to test diffusion in the 
models. 
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3. The fitting method 

The calibration of a binary system consists of adjusting the stel- 
lar modelling parameters so that the models reproduce the ob- 
servational data of the stars. The effective temperature and lu- 
minosity of a model, for a given mass and fixed input physics 
(EOS, opacity, nuclear energy generation rates), depend on the 
modelling parameters: the age f*, initial helium abundance 7 
and mixing length parameter a*. 

The effect that the variation of the initial helium abundance 
Yj and mixing length parameter a* has on the evolutionary 
tracks can be seen in Fig.^ An increase in a* moves the evolu- 
tionary tracks in the HR diagram towards larger effective tem- 
peratures at approximately constant luminosity, for stars with 
convective envelopes. An increase of the efficiency of the con- 
vection deepens the base of the convective envelope, while it 
does not affect the core. The temperature of the base of the ex- 
ternal convective zone therefore increases, and the adjustment 
of the adiabat to a higher temperature at its base causes the ef- 
fective temperature to increase, leaving the luminosity almost 
unchanged. 

Fig. |2] shows the evolution of the convective zones of 
EK Cep A and B. The convective envelope of EK Cep A disap- 
pears as the star approaches the main sequence, hence the small 
effect of the variation of a+ in the tracks of EK Cep A on the 
later phases of PMS evolution. 

A larger helium abundance shifts the evolutionary tracks 
on the HR diagram towards higher effective temperature and 
higher luminosity (see lower panel of Fig.Q. That is because a 
larger helium abundance causes the mean molecular weight to 
increase, which decreases the pressure. To maintain quasi-static 
equilibrium, the temperature must increase, which causes an in- 
crease in the energy generation rate and, thus, an increase in lu- 
minosity. This increase is more obvious in main sequence stars 
because of the large power-law dependence of the nuclear reac- 
tion rates on temperature. On the other hand, the gravitational 
energy generation rate increases only linearly with temperature 
and the rise in luminosity due to the larger helium abundance 
is smaller for PMS stars. 

We can, then, produce a stellar model with a higher ef- 
fective temperature by increasing the parameter a*; to change 
only the luminosity of the model, we can change simultane- 
ously Y\ and a*. To obtain observables as close as possible to 
the observations, we use the x 2 fitting developed by Lastennet 
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Fig. 1. Upper panel: variation of the evolutionary tracks with 
the mixing length parameter a. Full lines are evolutionary 
tracks for or=1.3; dashed lines are for a=l.6; and dotted lines 
are for a-\ .9 Lower panel: variation of the evolutionary tracks 
with the initial helium abundance. Full lines are evolutionary 
tracks for 7=0.24; dashed lines are for 7=0.26; and dotted lines 
are for 7=0.28. 



et al. ( 1999 see also Morel et al. 2000b I. It corresponds to min- 
imising the following functional: 



X 



Z 



log T m ° d 

b eft,* 



log r, 



eff,* 



Cr(log r ett> ) 



j ogZ mod_j ogU 



cr(log L*) 



(1) 



This functional depends on the modelling parameters through 
log T™^ and logL™ od . For a grid of modelling parameters we 
have computed the evolution of models with the masses and 
metallicities of EK Cep A and B. The functional^ 2 was com- 
puted for each of these models using the above equation. The 
parameters that minimised^ 2 were selected as the solution. 

4. Results 

The modelling parameters that minimise \ 2 are given in 
Table|2j together with the same parameters for a calibrated so- 
lar model using the same input physics. The uncertainties are 
derived by changing one parameter at a time around the cali- 
brated value, keeping all the others constant. The value for qta 
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Table 2. Modelling parameters for a model of EK Cep ly- 
ing within the uncertainty boxes for logT^* and logL+/LQ 
(see text). The parameters for the Sun are for a calibrated so- 
lar model using the same input physics. We do not present the 
value for or a because there is little dependence of the location 
of the model of EK Cep A on this parameter (see Fig.[Q. 



EK Cep 



Sun 



a B 1.363*$; 


1.632 


Y 261+ 0005 


0.265 


U 26.8!}| Myrs 


4.52 Gyrs 



G 

*S, ^ 7 EK Cep A 
>-3 



2, b 



3.9 



EK Cep B 



3.8 
lQ glO T e,f 



3.7 



3.6 




t, (Myr) 



t. (Myr) 



Fig. 2. Convective zones and radii of EK Cep A (left panel) and 
B (right panel). Vertical lines represent the extent of convective 
zones; the vertical dotted line in both panels represents the best 
model, with an age of 26.8 Myrs. 
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Fig. 3. Evolutionary tracks of EK Cep A and B on the HR di- 
agram. The uncertainty boxes of both stars are also shown. In 
the lower left panel, an enlargement of the locus of EK Cep A 
is shown; the arrow shows the location on the HR diagram of 
the calibrated model. 



is not listed in Table |2] since the dependence of the evolution- 
ary tracks on this parameter is very weak (see section[3] and in 
particular Fig.Q. We used a a = 1.340 to produce these mod- 
els. 

Our calibrated value of a# is lower than the value of a* 
for calibrated models of solar-mass stars (see, for example, 
Morel et al. 2000b for the case of a Cen, as well as the value 
for the solar case we present in Table Typically, one has 
a* -1.6-1.8 for calibrated solar-type stars. D'Antona & 
Montalban (2003 1 argue that in order to reproduce both the 
lithium depletion and the location of PMS tracks in the HR di- 
agram, convection must be highly inefficient during the PMS; 
that is, it must have a very low a* ~ 1. Our value of ckb is 
between that appropriate for a calibrated MS model and the 
value argued by D'Antona & Montalban for PMS models when 
they are depleting their lithium. This could reflect the fact that 
EK Cep B is in an evolutionary phase between the main lithium 
depletion phase (which occurs for stars with this mass around 
the time a radiative core starts to develop) and the ZAMS. 

The sequences of models for the evolution determined with 
the best set of parameters are shown in Fig.|2]and Fig. [3] Both 
components of EK Cep have a convective core; EK Cep A is 
a main sequence star with an energy production mostly due to 
the CNO cycle, while EK Cep B has just acquired a convective 
core due to the start of 12 C burning. 

This convective core in EK Cep B will be lost shortly after 
the arrival on the ZAMS. The ZAMS occurs at an evolutionary 



age of ? Z ams,A-10.8 Myrs and f za ms,B=51.3 Myrs for EK Cep A 
and B, respectively. 

From our estimated age of around 27 Myrs, follows that 
EK Cep A is already in the main-sequence, while EK Cep B is 
indeed a PMS star (see Figs|2]and[3}. The fact that both com- 
ponents are in different stages of evolution allows for a precise 
age determination, as the solution minimising Eq.^is strongly 
dependent on 

5. The problem of the radius 

Our best models for EK Cep A and B do not fit, however, the 
observed radii. Since EK Cep is an eclipsing binary, the radii 
of both stars are known with great precision; the fact that these 
models do not fit the radii is, therefore, a serious problem. 

We can fit the radius of EK Cep B to the observations using 
a lower in the models. Fig. 0] shows the evolutionary tracks 
for EK Cep A and B on the logR/RQ-logL/LQ plane, using 
as = 1 .28 (all other parameters being the same as in the previ- 
ous section). However, the same cannot be done for EK Cep A 
since the model depends weakly on the mixing length parame- 
ter a. 

We performed the fitting decribed in section[3]using the ef- 
fective temperature and radius as the observables. The param- 
eters that minimise^ 2 are now Y[ = 0.331, f* = 13.9 Myrs 
and as = 168. The helium abundance seems very large for the 
metallicity of EK Cep, but it reproduces the results presented 
in Lastennet et al. (2003 1 for the binary UV Psc. However, a 
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Fig. 4. Evolutionary tracks of EK Cep A and B on the 
logR/RQ-logL/LQ plane. These tracks were calculated using 
as = 1.28. Crosses show the location of models with 26 Myrs. 

Table 3. Characteristics of the best models, with Y\ = 0.261, 
as = 1 .28 and f* = 26 Myrs. T c and p c are the central densities 
and temperatures, given in 10 6 K and g cm 4 , respectively. R co 
is the radius of the convective core and R ce the radius of the 
base of the convective envelope. The same quantities are also 
given for a calibrated solar model using the same input physics 
(L Q = 3.846xl0 33 ergs s _1 and R Q = 6.9599xl0 10 cm). 





EK Cep A 


EK Cep B 


Sun 




21.11 


12.81 


15.45 


Pc 


73.96 


77.13 


149 




8974 K 


5550 K 


5777 K 


ljlq 


15.74 


1.465 


1.000 


RJRq 


1.644 


1.311 


1.000 




0.120 


0.0536 




Kc/R* 


0.993 


0.751 


0.729 



model of EK Cep A constructed using these parameters still 
has too large a radius. 

Fig- ID shows the evolutionary tracks in the HR diagram 
of models of both components calculated using the parame- 
ters given in the previous paragraph. Since these new param- 
eters fail to improve considerably the fitting to the radius of 
EK Cep A, while the location in the HR diagram of the model 
of this star falls well outside the error boxes, we prefer the solu- 
tion given in Table[2]except for ob ', to fit the radius of EK Cep B 
within the observational uncertainties, we choose = 1.28. 
The age that best reproduces the observations is f* = 26 Myrs. 
The characteristics of the best models, calculated using these 
parameters, are shown in Table[3] 

The fact that models of EK Cep A depend so weakly on 
the mixing length parameter a corresponds to having, in fact, 
only three "effective" modelling parameters, Y{, a% and f+. 
Thus, to fit the observables we need an extra parameter, related 
to the description of the interior physics affecting mainly the 
high mass component. This could be the differencial rotation 
of EK Cep A, as described in Yildiz ( 2003 1. 
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Fig. 5. Evolutionary tracks of EK Cep A and B on the HR 
diagram. These tracks were calculated using Y[ = 0.331 and 
org = 1 .68, and represent the best fit to the radii of both stars. 
Crosses show the location of models with 13.9 Myrs. 

6. The role of the time step 

Given the importance of EK Cep B in the age determination, it 
is fundamental that its evolutionary phase (PMS) is adequately 
calculated. Consequently, we have performed a detailed analy- 
sis on the prescription that should be adopted for defining the 
integration time step in the PMS evolution of this binary. 



EK Cep A 




EK Cep B 




4 3.98 3.96 3.94 3.92 
log,. T.„ 



3.78 3.76 3.74 3.72 
log,. T „ 



Fig. 6. The role of the time step. Left panel: EK Cep A; 
in full, the evolutionary track with a maximum time step of 
0.025 Myrs; dashed, the evolutionary track calculated using the 
time step prescription appropriate for the MS evolution. Right 
panel: the same for EK Cep B. In full, the evolutionary track 
with a maximum time step of 0.1 Myrs. 

Figure [6] shows again the evolutionary tracks in the HR di- 
agram for both components of the EK Cep system. The effects 
of an inappropriate choice of the time step are shown; these ef- 
fects are only relevant in the PMS phase because in this phase 
the star derives its energy from gravitational contraction; the 
rate of gravitational contraction is incorrect if the time step 
between two consecutive models is too big. The evolutionary 
tracks represented by a dashed line in Fig. |6] were calculated 
using a time step prescription that is appropriate for the MS 
evolution: the time step for the next model is estimated based 
on the evolution of the chemical composition. That is, the faster 
the evolution of the chemical composition, the shorter the time 
step used. This prescription is not useful for models in the PMS 
phase, since for most of the evolution the chemical composition 
does not change considerably. Instead, we must use a prescrip- 
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tion that does not allow the structure of a star to change too 
much between consecutive time steps. 

We can gauge how much the star has contracted between 
two consecutive time steps by calculating the gravitational en- 
ergy generated between them. To do this we must first calculate 
the model corresponding to f+Af. If the gravitational energy 
generated between the two time steps (E g ^-TAs, where s is 
the specific entropy and T is the temperature) is bigger than a 
given maximum value, a new model must be calculated corre- 
sponding to f+Af for a smaller Af. There are, however, some 
numerical difficulties in this approach if the maximum gravita- 
tional energy allowed to be generated between two consecutive 
time steps is too small. 

We have chosen a different prescription instead, by using 
the time the star spends in the PMS phase (fpMs)- This time 
is expected to scale, roughly, with the Kelvin-Helmholtz time 
scale, fpMS-fKH-GM^./2/?*L+, where and L+ are the radius 
and luminosity of the star on the main sequence. These values 
can be replaced by their ZAMS values since for most of the 
time the stars are in the radiative part of the PMS evolutionary 
tracks where and L* do not change very much. 

For a given mass range, R+ocM^ and L+ccMl, on the main 
sequence. In the mass range of interest here, ^0.6 and /p3. 9 
(see Kippenhahn & Weigert 1991 ch. 22), it follows that 

fPMS-f^-M^-^M, 2 - 5 . (2) 

We scale the maximum time step allowed with the time the star 
spends in the PMS phase, so that the evolutionary tracks in the 
PMS phase are equally divided. For a IMq star, we have con- 
sidered different values for the maximum time step; for values 
shorter than about 0. 1 Myrs, the evolution of the star does not 
change for different values of the time step. Accordingly, the 
maximum time step we allow in the evolution during the PMS 
phase is considered to be given by 

(.» v — 2,5 
S§) Myl ' S ' (3) 

(see also Kippenhahn & Weigert 1991 ch. 28). 

The effects of different prescriptions for the time step on the 
calibration of EK Cep can be seen in Fig. [6] The locus of the 
calibrated model of EK Cep A (left panel) in the HR diagram 
does not change if we do not control the time step carefully 
because the star derives most of its energy from nuclear reac- 
tions. On the other hand, the locus of the calibrated model of 
EK Cep B changes considerably; the luminosity of the model 
calculated without the restriction on the time step is apprecia- 
bly lower. In fact, we can not fit the observations within their 
error boxes for both stars if we do not restrict the time step 
in this way; the luminosity (and radius) of the low mass com- 
ponent is too low. This partially explains the results cited in 
Martin & Rebolo ( 1993 1, that the luminosity of EK Cep B was 
slightly higher than predicted by pre-main-sequence models. 
The radius of our model of EK Cep B (with a large time step) 
is below the error box in radius, a result also obtained by Claret 
et al. (1995}. 

The evolution of the models is faster if we do not restrict 
the time step. Consequently, instead of the value we found 



(26.8 Myrs) the calibrated models of EK Cep would just have 
an age of =2 23 Myrs when using the same prescription as for 
the MS evolution. 

7. The effect of overshooting 

We also calculated evolutionary models including the effects of 
convective overshooting. Convective overshooting is due to the 
inertia of the convective elements; since they have a non-zero 
velocity when they reach the border of the convective zone, as 
given by the instability criteria, they overshoot it. When they 
reach a radiative zone, they become colder than the surround- 
ing material (since in a radiative zone V < V ac |) and are there- 
fore heated by their surroundings: they transport energy down- 
wards. To transport all the energy upwards, the energy trans- 
ported by the radiation must be greater than the total luminos- 
ity, that is, V > V la d. The effect of the overshooting is, then, 
to fix the temperature gradient equal to the adiabatic gradient 
and to extend the mixed zone by a certain amount. We used the 
following prescription for the extent d m of the overshooting in 
the core, 

d ov = a m x Min (H P , r co ) , (4) 

where a ov is a free parameter and r co is the radius of the con- 
vective core. 




Fig. 7. Evolutionary tracks of EK Cep A and B on the HR di- 
agram. These tracks were calculated including the effects of 
convective overshooting, with ff ov = 0.20. 

Fig. shows evolutionary tracks calculated using a m = 
0.20. The effects of the overshooting on the tracks are small 
for the track corresponding to EK Cep B; it merely extends the 
PMS phase, and makes the drop in luminosity near the end of 
the PMS phase more abrupt. This is caused by the fact that the 
central convective zone that appears at an age of about 25 Myrs 
(see again Fig.fJJ has a higher amount of 12 C to burn (the mixed 
zone extended by overshooting has a higher mass); the energy 
produced by 12 C-burning is higher, which causes the luminos- 
ity to drop more abruptly, and the fact that there is more I2 C to 
burn means that the ZAMS will arrive later (by about 10 Myrs; 
the precise value depends on the parameter or ov ). 
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On the other hand, the track corresponding to EK Cep A 
is considerably changed with the inclusion of overshooting; it 
displays an extra "loop" after the first drop in luminosity (com- 
pare Figs.[3]and0i. This "loop" in the evolutionary track in the 
HR diagram is also seen (but not discussed) in the evolutionary 
tracks of Siess et al. ( I2000> . It is caused by the following: as 
the 12 C in the central convective zone is been burned, its abun- 
dance drops and so does the rate of burning. As can be seen 
in Fig. 13 the central convective zone recedes towards the cen- 
tre, leaving behind a zone of partially depleted 12 C. When the 
burning of I4 N becomes more efficient due to the rising central 
temperature, the convective zone expands again, englobing in 
the mixed central zone the unburned 12 C at the external regions. 
Now, as will be discussed in a subsequent paper (Marques et al. 
in preparation), the extent of the mixed zone due to overshoot- 
ing leaves more 12 C unburned in these regions; when this sec- 
ond expansion of the convective zone engulfs fresh 12 C, it does 
so at a rate faster than the burning time. The rising of the abun- 
dance of 12 C in the central convective region (which has now 
a much higher temperature than it had in the previous 12 C- 
burning phase) causes a large increase in the energy produc- 
tion through nuclear reactions in the core, which in turn makes 
the central convective region increase very rapidly, with the 
consequent rapid drop in the luminosity of the star. This ex- 
tra "loop" in the evolutionary track in the HR diagram takes 
some 1.5 Myrs; thus, the arrival on the ZAMS is delayed by a 
significant amount (for a m = 0.2, about 2.2 Myrs). 

For the present work, all these effects are of little impor- 
tance. The age of the system EK Cep is greater than the age 
at which these effects take place for a star with the mass of 
EK Cep A, while it is lower than that at which overshooting 
starts to make a (small) difference for stars like EK Cep B. Our 
calibration is not, therefore, changed by the inclusion of con- 
vective overshooting. An alternative for studying the presence 
of overshooting in this system is the use of asteroseismology 
(e.g. Monteiro et al. l2000> . 

8. Conclusions 

EK Cep is an excellent candidate to test stellar evolutionary 
models because it has one member in the MS phase and the 
other in the PMS phase. This difference in the evolutionary 
phase of the components results in a better calibration of the 
binary age of 26.8 Myrs. 

We have shown that an incorrect treatment of the time step 
for the evolution yields models with a lower luminosity during 
the PMS phase. If both stars are in the PMS phase, a large time 
step can be compensated for by a higher initial helium abun- 
dance (to increase the luminosity of the models; see Fig. Q. 
This cannot be done if the stars are in different evolutionary 
phases because the increase in the initial helium abundance will 
make both stars brighter, while only the PMS component has a 
lower luminosity, due to the larger time step. The dependence 
of the calibration on this key aspect of the modelling might 
explain the results obtained by Claret et al. ( 1995 1 and other 
results cited in Martin & Rebolo ( 1993 1, in particular that the 
radius (and luminosity) of EK Cep B are underestimated by 
evolutionary models. However, we cannot reproduce the radius 



of EK Cep A with our models. The presence of a fast rotating 
core (as in Yildiz 2003 (would improve the situation. The use of 
a correct time step solves only part of the problem, by increas- 
ing the radius of EK Cep B and reducing the need for a larger 
radius of EK Cep A. We stress that our models were calculated 
using the metallicities obtained by Martin & Rebolo (1993); 
if a new determination of the metallicity is made, it would be 
possible to discern more clearly the origin of this inconsistency. 

In the light of the new calibration, it becomes clear that 
EK Cep is an excellent binary to test different aspects of the 
physical processes of the evolution. Given the particular con- 
figuration of this system, one component (the primary) is insen- 
sitive to the mixing length parameter and age, while providing 
a very precise constraint on the helium abundance for the sys- 
tem. The effect of the overshooting on the evolutionary tracks 
of EK Cep A, although not relevant for this work, is important 
near the end of the PMS phase. We will explore the effects of 
the overshooting during the PMS evolution for various types 
of stars in a subsequent paper (Marques et al. in preparation). 
EK Cep B, still in the PMS, provides a very accurate indica- 
tion of the age of the system, while being independent of the 
overshooting prescription. Such a combination makes EK Cep 
a test case of how the detailed observation of young binaries is 
of great relevance for improving the modelling of the evolution 
in this phase. 
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